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Abstract 

At LEP II it is hoped to measure the W mass to an accuracy 
of around 40 MeV. This will require direct reconstruction of the 
mass of the W from its decay products in both the semi-leptonic 
and hadronic decay channels. Final state perturbative reconnec- 
tion effects in hadronic decays are considered and their effect on 
6-jet distributions and the reconstructed mass. The perturbative 
mass shift is found to be ~ 50 keV in the negative direction. 
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1 Introduction 



One of the main goals of LEP II will be an accurate determination of the 
mass of the W boson. An integrated luminosity of 500 pb™ 1 suggests that 
an accuracy of 30 — 50 MeV [|TJ could be reached. The process 

e + e~ — > W + W~ — ► 4 fermions (1) 

can be split into three distinct classes depending on the type of decay of each 
W-boson. 

• Purely leptonic. Both W-bosons decay to leptons. There are two neu- 
trinos and reconstruction of the event from observed charged lepton 
momenta is not possible. Branching ratio for this channel ~ |. 

• Semi-leptonic. One W-boson decays to leptons, the other decays hadron- 
ically. One neutrino is produced, but the missing momentum can be 
reconstructed using energy-momentum conservation and assumptions 
about the initial state radiation. Branching ratio for this channel ~ |. 

• Fully hadronic. Both W-bosons decay hadronically. All momenta are 
observable. The momenta directions are well resolved, while the en- 
ergy resolution can be improved via kinematic fits (ie imposing the 
constraints of energy and momentum conservation). Branching ratio 
for this channel ~ |. 

In order to achieve the greatest accuracy the W mass must be recon- 
structed using both the semi-leptonic and the fully hadronic decay channels. 
However Ws decay very rapidly so one expects that the space-time sepa- 
ration of the two decays should be ~ 0.1 fm. This is small compared to 
the typical scale of hadronization ~ 1 fm, thus in the case of fully hadronic 
decay there are two evolving hadronic systems with considerable space-time 
overlap. There is the possibility that the two systems do not evolve indepen- 
dently, but influence each other. 

These influences fall into two categories^ - Bose-Einstein correlations be- 
tween identical bosons in the final state (typically pions) 0, El M, and a 

2 I neglect the effects of electroweak interactions between the two systems as these have 
been considered elsewhere^], || 
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Figure 1: e + e — ► W + W —>■ qqQQ. The blob represents a sum over the 
three lowest order production amplitudes. 

re-arrangement of the colour flow of the evolving systems at either the per- 
turbative or hadronization level[|7|, ||, |9[ pL0f| - There has been much work 
on the effects of colour re- arrangement at the hadronization level, however 
hadronization is poorly understood and progress can only be made through 
constructing models. It is interesting to note that the models of colour re- 
connections in the hadronization phase give rather varied predictions |], |Tl| 
for the effects on physical observables such as mean charged multiplicity or 
reconstructed W mass, and so such measurements may probe directly aspects 
of the confinement mechanism. 

In this paper I will examine the effects of colour reconnection at its lowest 
non-trivial order in perturbation theory. In section two I will explain why 
these effects should be small and how they can be calculated directly. In 
section three I shall present results for the effects of colour reconnection on 
various distributions including the W mass. The conclusions will be found 
in section four. 

2 Perturbative Reconnection 

Perturbative reconnection appears as higher order corrections to the process 
shown in fig 1, in which gluons are exchanged between the evolving quark 
systems. A possible reconnection diagram is shown in fig 2; here a real 
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Figure 2: A possible interference term involving one exchanged gluon (the 
shaded blobs represent a sum over the three W-pair production amplitudes). 

gluon emitted from one decay system interferes with the similar emission 
from the other decay system. One may also consider the analogous virtual 
interference corresponding to the exchange of a virtual gluon between decay 
systems. Within perturbation theory these interference terms are zero due 
to colour conservation?], and so one must consider the exchange of at least 
two perturbative gluons. 

A full calculation of the 0{a 2 s ) corrections is beyond the scope of this 
paper, however it is possible to examine QCD interference effects in the 
production of 6 jets|f[2| via 

e + e~ — > W + W~ — > qqqqgg (2) 

in which interferences appear between the lowest order diagrams. Two pos- 
sible interference terms are shown in fig 3. 

These diagrams contain only two colour loops, compared with the di- 
agrams for gluon emissions within each decay system which contain four 

3 However it is not impossible for a colour octet to be exchanged between the two decay 
systems at the perturbative level, only to be balanced by a non-perturbative exchange 
in the hadronization phase. Such interplay between perturbative and non-perturbative 
connections is not considered here. 
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1 I 

Figure 3: Interference terms involving the exchange of two gluons (the W-pair 
production parts of the diagrams are omitted for clarity) 



loops. Therefore the interference terms are suppressed relative to the leading 
emission by -jkr where N c is the number of colours. 

There is further suppression due to the width of the W. Gluons radiated 
within a decay system are free to have any energy up to ~ M\y without 
pushing the W Breit-Wigner propagators off resonance. However gluons 
radiated between the decay systems (interference terms) must carry energy 
less than ~ IV or at least two of the W propagators must be pushed off 
resonance and that term will become suppressed. It has been shown that for 
inclusive quantities, where the I.R divergences cancel between real and virtual 
diagrams, that this leads to a suppression of perturbative reconnection effects 



by 0(-j^) [|13| , 15| . A rough estimate of the size of perturbative reconnection 



effects in W-pair production is thus: 

' U lO" 4 (3) 



Act a 2 s I\ 



and so one may estimate the possible mass shift as 



few MeV (4) 

N c 

This should be regarded as an order of magnitude estimate only. It is 
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Figure 4: Components used to calculate lowest order e + e~ — ► W + W~ — > 
qqqq. The shaded blob represents the sum of the three interfering amplitudes, 
and the WW propagators have been absorbed into the definition of T^ u 



clearly desirable to calculate experimental distributions in fixed order per- 
turbation theory and examine how they may be distorted by the effects of 
colour interference. 



Using the helicity methods of |L5] it is possible to construct the amplitudes 
for all the doubly resonant diagrams contributing to the qqqqgg final stateQ; 
there are 72 diagrams in total. In this method each amplitude is built up 
from relatively few component pieces that can be calculated separately. For 
example the amplitude for fig 1 can be written 



M (n, r 2 ; pi, p 2 , p 3 , p 4 ) = T^(r 1; r 2 ; pi +p 2 , p 3 +Pi)i^ ] (puM^ Os, Pa) (5) 

where the terms are defined in fig 4 and the momentum labels refer to fig 
1. The computational complexity is reduced by assuming massless electrons 
and massless quarks which has been done throughout this paper. 

In this way all of the amplitudes can be built up from just the production 
tensor T^ v and some 'decay currents' j™. In order to calculate all amplitudes 
for double gluon radiation four additional decay currents (fig 5) are needed. 
Note one must distinguish between and which are related by = 
jjp(ki <-> k 2 ) in order to obtain the correct colour factors. 

The decay currents may be contracted onto the production tensor to 
obtain eight distinct Lorentz-colour structures (colour-matrices are omitted 
for clarity) 

4 Strictly speaking this is not a gauge invariant set of diagrams, however one may see 
that a change of gauge leads to singly resonant contributions which are neglected. The 
amplitudes were evaluated in the unitary gauge. 
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Figure 5: Additional components needed to calculate the doubly resonant 
contributions to e + e~ — > qqqqgg. The shaded blob represents a sum over 
possible attachments of gluons to the spinor line which preserve the order of 
attachments. 
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Mi = T^jj 2 \p 1 ,p 2 ,k 1 )ji 2 \p 3 ,p 4 ,k 2 ) 

M 2 = M 1 (k 1 ^k 2 ) 

M 3a = TrjWfa^kuWPfapt) 

M 3b = T^ 3 f\p u p 2 MM)3 ( v\p^) 

M 4a = T^j^(p 1 ,p 2 )ji 4 \p 3 ,p 4 ,k 1 ,k 2 ) 

M, b = T^jW(p 1:P2 )jW(p 3 ,p, : k 1: k 2 ) 

M 5 = T^jf{p l ,p 2 MM)3 ( '\p^PA) 

M 6 = T^ti\ Pl ,p 2 )ti\p^p A MM) (6) 

then the total production amplitude is (suppressing colour matrices) 

M = Mi + M 2 + M 3a + M 3b + M 4a + M 4b + M 5 + M 6 (7) 

after squaring and summing over colours it is convenient to separate the 
squared matrix element into different parts depending on the form of the 
Breit-Wigner resonances. In this way one finds six distinct terms 



]T \M\ 2 = M 1 + M 2 + M 3 + M A + M 5 + M 6 (8) 

colours 



where 



M 1 = N%C%MiM{ (9) 
M 2 = M 1 (ki <-> k 2 ) (10) 
M 3 = N 2 C 2 F (M 3a M* 3a + M 3b M* 3b ) + N^C F M 5 M* 5 - 

l -N c C F {M 3a Ml b + M 3b M* 3a ) + 

l -N*C F {[M 3a Ml + M 5 M* 3a ] - [M 3b M* 5 + M 5 M* 3b ]) (11) 
M 4 = N 2 c CUM, a Ml a + M 4b Ml b ) + N 3 c C F M 5 M* 5 - 
l -N c C F {M ia Ml b + M 4b M* Aa ) + 

l -N*C F {[M, a Ml + M 5 M*J - [M 4b M* 5 + M 5 M* 4b }) (12) 
M 5 = -NcC F (MiMZ + M 2 Ml) (13) 



M 6 = ^N c C F [(M 3a + M 3b )(M Aa + M Ab )* + 

(M u + M ib ){M 3a + M 3b )*\ (14) 

In the above expression the terms Mi, M 2 , M 3 , M 4 denote the unrecon- 
nected parts, while M 5 and M 6 correspond to interference between the two 
decays. The separation into reconnected and unreconnected parts is (QCD) 
gauge invariant. 

One can find small regions of phase space where the interference contribu- 
tion to the total transition probability is as large as 10%, even for relatively 
energetic gluons (> 5 GeV). It is therefore not impossible that certain 6 jet 
distributions could be significantly distorted by perturbative reconnection. 

A Monte Carlo program was written to generate WW events with the 
qqqqgg final state using the Multichannel approach |TB| with 24 channels 



based on the kinematic properties of the contributing diagrams. Events were 
generated at a centre-of-mass energy of 192 GeV, although reconnection ef- 
fects are expected to be insensitive to the centre-of-mass energy within the 
LEP II range. In addition specific phase space parameterisations were com- 
puted which allowed efficient integration of both types of interference term. 

Six jet final states were defined according to a minimum invariant mass 
between partons. A lower limit of s cut = 1 GeV 2 was used. Strictly speaking 
this is too small for fixed order perturbation theory to be applicable, however 
the philosophy is that the results obtained will provide an upper limit on 
reconnection effects, since moving to larger values of s cu t generally reduces 
any effect. The six partons were clustered to four jets using the Durham 
algorithm [H]. Distributions for the Bengtsson-Zerwas angle [|18|j {xbz), the 
modified Nachtmann-Reiter angle fl9j (9nr) and the angle between the two 



lowest energy jets (034) were computed with and without the interference 
terms. These angles are defined by equation 15 below in which pi,P2,P3 
and P4 are the energy ordered jet 3-momenta. 



cos(x 



BZ 1 



(Pi X P2) ' (P3 X p 4 ) 
I Pi X P2| |P3 X P4| 



la \ IPl ~ P2j ' IP3 - P4j 

cos(9 NR ) = — r — (15) 

Pi - P2IIP3 - P4l 

With four jets there are three ways of pairing them. For each pairing 
the average of the invariant masses was computed. Thus for each event one 
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has three mass values corresponding to each of the three possible pairings. 
The mass closest to the input W mass was chosen as the mass estimate 
(this is only one of several possibilities suggested in 0). Distributions for 
the mass calculated in this way were also produced with and without colour 
interference terms. 

The difference between distributions with and without the interference 
terms was computed. This shows the distortion induced by the interferences. 

Finally the integral of the absolute value of the interferences was found 
for gluon energies greater than 2 GeV, 5 GeV and 10 GeV. These quantities 
are finite since the interference terms contain no collinear singularities (apart 
from integrable ones when three partons become collinear), and provide an 
indication of the possible size of interference effects in events with jet energies 
greater than 2, 5 and 10 GeV. 

3 Results 

The mean mass can be calculated with and without the reconnected terms. 
The result one finds depends on the choice of invariant mass cut, but must 
tend to zero as s cu t — > since the unreconnected terms are more singular than 
the reconnected terms in this limit. Mass shifts for a variety of invariant mass 
cuts on the final state are shown in the table below. 



Scut/GeV 2 


0.1 


1.0 


10.0 


100.0 


5M w /MeV 


-0.030 


-0.045 


-0.025 


~ -0.015 



The exact numbers are also slightly dependent on the reconstruction scheme 
used for defining the experimental W mass. 

Figure 6 shows the distribution of reconstructed mass using only the 
unreconnected parts of the matrix element (solid line). The dashed line 
shows one thousand times the change induced when the reconnected terms 
are present. The distributions for the mass under the full matrix element 
and unreconnected terms only differ essentially by a multiplicative constant 
of order 1.001. The mean value of the mass distribution is shifted by less 
than a part per million due to the presence of reconnected terms. 

Figures 7,8 and 9 show similar plots for the distribution of the Bengtsson- 
Zerwas angle, the angle between the two lowest energy jets a 34 and the 
Nachtmann-Reiter angle. It will be seen that the effect is at or below the 
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Mass/GeV 



Figure 6: Mass distribution for unreconnected events (solid line), change 
induced by reconnected terms xlOOO (dashed line), in arbitrary units 
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0.5 1 1.5 2 2.5 

Bengtsson-Zerwas Angle 



Figure 7: Distribution of Xbz, the Bengtsson-Zerwas Angle for unreconnected 
events (solid line), change induced by reconnected terms xlOOO (dashed line) 



per mille level and is essentially just multiplicative, distortions of the distri- 
butions occur at a much lower level. These effects can be understood within 
the soft interference limit. In the soft limit one may describe gluon radiation 
using eikonal vertices and the matrix element squared becomes. 

\M\ 2 = \M \ 2 (H(h, k 2 ) + AG(h)G(k 2 )) (16) 

where H(ki,k 2 ) is the soft unreconnected distribution, A is some constant 
that will depend on the energy resolution and W width. G{k) is the re- 
connected distribution (note that at this order the reconnected gluons are 
radiated independently, however this is not true in higher orders) and is 
given by 



G(k) = ' P ^ + ^ 2 " Pz) (pi • Pa) (gg • Pa) 

(pi • k) (p 4 • k) (p2 • k) (pa • k) (pi • k) (p 3 • k) (p2 • k) (p 4 • fc) 

(17) 

One may integrate over the directions of each emission to find the en- 
hancement due to soft interference between decays: 
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Cos(.alpha34) 



Figure 8: Distribution of 005(0:34) for unreconnected events (solid line), 
change induced by reconnected terms x500 (dashed line), in arbitrary units 



ItiJ-W^ NH (18) 

\(Pi ■P3){P2 -Pa)) 

where the momenta are as defined in fig 1. 

The effect of the reconnection terms is essentially to enhance coplanar 
configurations where some invariant masses can be much larger than others. 
In configurations where the W decay planes are at right angles, none of the 
parton directions can become close and so the argument of the logarithm 
in equation (18) is close to one and there is little enhancement. In most 
approximately coplanar configurations the BZ angle will be close to either 
or 7T as both p x x p 2 and p 3 x p 4 (energy ordered momenta) are likely 
to point out of the decay plane and hence be either parallel or anti-parallel. 
Thus one expects enhancement around these values. 

The situation for 003(034) is not quite so straight forward. A similar 
argument favours 003(0:34) ~ 1 however this configuration is suppressed by 
the jet reconstruction kinematics; one would need two low energy quarks 
and both gluons radiated in approximately the same direction and to be 
clustered as two distinct jets. However the configurations corresponding to 
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Cos(NR) 



Figure 9: Distribution of cos(#7v.r), the Nachtmann-Reiter Angle for unre- 
connected events (solid line), change induced by reconnected terms x500 
(dashed line) 
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Figure 10: Kinematic configurations which are a) Enhanced by interference 
but kinematically suppressed (cos(a34) ~ 1), and b) Enhanced by interfer- 
ence (cos(a 3 4) ~ —1), since (pi-p4)(p2-P3) > {Pi'Pzjipi'PA}- The momentum 
labels refer to fig 1, and in this example p 2 and p 3 correspond to the lowest 
energy jets. 

cos(o;34) ~ — 1 can be enhanced (see fig 10). 

A similar argument for 9nr is not so apparent as its geometrical interpre- 
tation is less clear (the angle between the axis defined by the vector between 
the two lowest energy jets and that between the two highest energy jets). One 
may construct the enhancement due to equation (18) and find qualitatively 
the same shape as observed in figure 9. 

The absolute value of the interference terms was integrated over the region 
defined by uj > 2 GeV, uj > 5 GeV and uj > 10 GeV where uj is the minimum 
gluon energy. This was done for s cut = 10, 1.0, 0.1, 0.01 GeV 2 and illustrates 
collinear finiteness (see results below). 



0"|mt|/pb 


Scut/GeV 2 


10 


1.0 


0.1 


0.01 


uj/GeV 


2 


0.015 


0.029 


0.032 


0.033 


5 


0.0057 


0.0082 


0.009 


0.01 


10 


0.0016 


0.0023 


0.0026 


0.0026 



The errors on these numbers are around 4% each. Note that these are the 
integrated absolute value of the interference terms, the actual contribution of 
the interference terms to the cross-section is typically an order of magnitude 
smaller due to large cancellations. 
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4 Conclusions 



Effects of perturbative reconnection are not necessarily small, however the 
regions of phase space in which sizable effects can occur are small. Most 
experimentally interesting distributions are unaffected by reconnection at 
the perturbative level apart from a multiplicative factor close to unity. In 
particular the mass distribution is shifted by less than one part per million 
by lowest order reconnection effects in 6-jet events. Distributions sensitive 
to soft momenta seem to show greater distortion, however these effects are 
well below the per mille level and so unlikely to be seen at LEP II. 

The integration of the absolute value of the reconnection terms for gluon 
energies above 5 GeV shows that the maximum effect could only be equiva- 
lent to a few events at LEP II and can probably be neglected at this level of 
statistics. 

Of course reconnection effects summed over higher terms, or within the 
hadronization phase need not be negligible and these effects still need to be 
addressed. 
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